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Abstract 

Half of the world's forest is in boreal and sub-boreal ecozones, containing large carbon stores and fluxes. Carbon lost from 
headwater streams in these forests is underestimated. We apply a simple stable carbon isotope idea for quantifying the CO2 
loss from these small streams; it is based only on in-stream samples and integrates over a significant distance upstream. We 
demonstrate that conventional methods of determining CO2 loss from streams necessarily underestimate the CO2 loss with 
results from two catchments. Dissolved carbon export from headwater catchments is similar to CO2 loss from stream 
surfaces. Most of the CO2 originating in high CO2 groundwaters has been lost before typical in-stream sampling occurs. In 
the Harp Lake catchment in Canada, headwater streams account for 10% of catchment net CO2 uptake. In the Krycklan 
catchment in Sweden, this more than doubles the CO2 loss from the catchment. Thus, even when corrected for aquatic CO2 
loss measured by conventional methods, boreal and sub-boreal forest carbon budgets currently overestimate carbon 
sequestration on the landscape. 
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Introduction 

Boreal and sub-boreal ecozones are large (about 1 7 million km^) 
and about half the world's forest cover [1]. Forests in these zones 
contain large carbon stores and contribute significant fluxes in the 
global carbon budget [2]. Part of the carbon fixed from the 
atmosphere by forests is returned to the atmosphere via aquatic 
surfaces [3-6]. In the past, this flux has been ignored in the 
construction of carbon budgets for forested watersheds and is now 
predicted to be higher than expected [7]. In Sweden, annual 
estimates of carbon sequestration in forests are 10% lower if 
aquatic carbon losses are included [8]; similar amounts are 
expected for Norway [9]. In the USA, stream and river CO2 loss 
may be five times greater than previously thought [10]. The 
difficulty in accounting for the aquatic loss lies in inadequate 
information on the enormous number of small headwater streams 
with high CO2 concentrations. In boreal Sweden, 90% of total 
stream length is in catchments less than 15 km^ [1 1]. 

Neglecting the aquatic export of terrestrially fixed carbon (via 
dissolved inorganic carbon (DIG), CH4, particulate organic carbon 
(POC), and dissolved organic carbon (DOC)) in landscape carbon 
budgets results in overestimating net ecosystem exchange by 25% 
to 70% [12-14]. The active pipe model for inland waters [15] 
increased estimates of total terrestrial export by two-fold while 
concluding that CO2 loss to the atmosphere was necessarily 
underestimated because small streams were entirely excluded for 
lack of emission and distribution data [3-7,10]. Subsequently, 'in- 
stream heterotrophy' was added to the stream portion of the 
budget [5] but other DIG sources, e.g., DIC-rich groundwaters, 
are poorly known and therefore excluded. The knowledge gap 



between the significant contribution of headwater streams in a 
basin and the small amount of data about them has been termed 
aqua incognita [11]. 

AU streams in an extensively studied 4th order Swedish boreal 
catchment were supersaturated in CO2 [16]. There, CO2 loss 
from stream surfaces directly to the atmosphere was 30%-50% of 
the dissolved carbon (DIC-I-DOC) exported from the forest 
(2.9 gC/mVyr of 9.8 gC/mVyr) [17,18]. In an Alaskan head- 
water stream network, greatest variability and mean fluxes of CO2 
were in the first order streams [19] with 9.0 gC/m^/yr lost from 
the entire Yukon River basin [20]. 

In sub-boreal catchments in Canada [21] and a peat catchment 
in northeast Scotland [13], CO2 loss from stream surfaces was 
36% and 34% of the dissolved carbon exported from the 
catchments. These studies noted that gas loss from stream surfaces 
to the atmosphere is significant, but all neglected the higher CO2 
evasion rates from stream surfaces upstream of their sampling 
locations. Without good estimates of CO2 fluxes, ecosystem-scale 
metabolism calculations also contain this uncertainty and bias. 
Die in headwater streams is a net result of a number of processes, 
including: dissolution of carbonate and weathering of some 
minerals in soils and bedrock, in-stream biotic respiration and 
fixation, exchange with atmospheric COg, and shallow ground- 
water input [17,22-25]. Although it is commonly assumed that 
CO2 saturation is controlled by mineral weathering and in-stream 
biotic activity in carbonate systems, boreal forests are, in many 
cases, characterized by high-DIC shallow groundwaters derived 
from terrestrially respired organic carbon [8,24]. 
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Here, we focus on silicate bedrock catchments because they 
constitute 45% of boreal catchments [1] and we avoid the 
complications of the contribution of significant amounts of DIG 
from carbonate [26]. Silicate-dominated boreal landscapes con- 
tribute significandy to global carbon fluxes [27], Headwater 
streams in these catchments have low to moderate pH values and 
moderate to high concentrations of CO2 that is modem in 
radiocarbon terms [28,29] . Thus, DIG from shallow groundwater 
is directly related to soil and root respiration instead of long-term 
mineral weathering. 

In undisturbed, shaded, nutrient-poor, sihcate-bedrock, boreal 
headwater streams, shallow groundwater is the major source of 
DIG since in-stream respiration is low relative to gas exchange 
with the atmosphere [30-33]. Furthermore, in sUicate terrain, 
shallow groundwater retains the <5''G value of the DIG input from 
DOG and soU OC decomposition and exhibits modern radiocar- 
bon ages [28,29] indicating weathering is unimportant at this 
scale. Thus for small streams of up to a few hundred metres length 
in this landscape, soil and root respiration is the key (5''G-DIG 
source. 

Degassing of GO2 from these streams causes the 5''G-DIG to 
increase due to well known equilibrium isotopic fractionation 
between GO2, HCO3", and G03^ and during GO2 loss to the 
atmosphere — since the (5'^G- CO2 value is less than the bulk 
^'^G-DIG value, the loss of GO2 wUl cause the remaining ^'^G- 
DIG value to increase [34,35]. Additionally, since GO2 loss does 
not affect carbonate alkalinity, there is a concomitant increase in 
pH as degassing occurs along the stream [36]. This CO2 loss also 
includes GO2 lost as groundwater transits the riparian zone 
adjacent to the stream since diese DIG and ^'^C-DIC changes are 
indistinguishable from those in the stream proper. Recent work 
presented a moderately complex modelled formulation of th(- 
degassing effects on ^'^G-DIC and focused on several French 
streams [35]. 

At a much broader level, we can use the characteristic degassing 
trajectories of decreasing DIG vs increasing ^'^C-DIC to estimate 
(1) how far along the degassing curve a particular sampling site in a 
stream is and (2) how much CO2 loss there must have been 
upstream of the sampling site. This gives us the opportunity to use 
discrete samples from a stream to estimate the upstream CO2 loss 
without direct knowledge of the DIG concentration of the 
groundwater end-member. Since groundwater (5'^G-DIG can be 
easily constrained based on G3 plants and soil organic matter 
formation [22,28,37] in silicate-dominated catchments, the aver- 
age groundwater DIG concentration can then be determined as a 
result of the modelling process. Discrete samples are available 
from the mouths of streams where whole-lake mass balances have 
been completed [38], along stream networks [16,39], and 
potentially from archived samples and data sets [40]. Because 
only in-stream samples are needed, this technique makes it 
possible to cover large catchment areas. Logistically, headwater 
stream CO2 fluxes can be rapidly assessed without the large 
investment required to study only one catchment in detail. 

Given the (1) poor knowledge of the extent of headwater stream 
surface area (lengths and widths [11,24], but see recent methods 
using digital elevation models [18,20] and branching theory [41]) 
and (2) narrow constraint on ^'^C-DIC in groundwater, we 
suggest these small streams can be studied with archived and new 
samples. We place the results in a context of landscape-scale 
degassing trajectories where (5'^C-DIC of groundwater, stream, 
and small lakes all differ. We focus on improving CO2 loss 
estimates from headwater stream surfaces in two catchments: (1) 



Harp Lake, Ganada where stream GO2 loss is situated in a 
catchment G budget and (2) Krycklan, Sweden where stream GO2 
losses are higher than previously estimated. 

Materials and Methods 

Site description 

Two field sites of contrasting vegetation and location were 
selected: Harp Lake catchment, Canada (45°22' N, 79°08' W) and 
Krycklan catchment, Sweden (64° 14' N, 19°50' E). Three 
headwater streams in the Harp Lake catchment were studied. 
AU are underlain by the gneissic bedrock of the Ganadian Shield 
[42] with mixed forests of sugar maple {Acer saccharum), beech 
{Fagus grandifolia), yellow birch (Betula alleghaniensis), white 
pine (Pinus strobus), and aspen (Populous tremuloides). Harp Lake 
catchment receives 1033 mm/yr precipitation. Mean annual air 
temperature is 5°G. Stream catchments ranged from 3.7 km^ to 
21.7 km^ with stream lengths from 170 m to 760 m. AU streams 
were first or second order and catchments had varying amount of 
wetland cover (<0.3% to 8.5%) [43,44]. 

Stream samples from around Harp Lake were collected from 
July 1989 to November 1991 across various flow conditions near 
the mouth of each stream (« = 1 1 1). DIG samples were collected in 
Pyrex culture tubes with polycone caps and analysed colorimet- 
rically as per standard methods in the Ontario Ministry of the 
Environment lab at Dorset, Ontario, Ganada. Samples for (5'^G- 
DIG analysis were collected in 500 mL glass bottles with polycone 
caps and analysed by standard off-line dual-inlet IRMS tech- 
niques. Groundwater was sampled from piezometers near the 
middle of the Harp 4-21 catchment (w = 8) [28,29]. An automated 
titrator (PG-Titrate Man-Tech Associates) was used to measure 
pH. Stream temperature was measured in the field. Stream DOC 
concentrations varied with stream flow with a median concentra- 
tion of 5.1 mgC/L (range of 1.1 mgC/L and 17.2 mgC/L) and 
flow- weighted annual average around 10 mgC/L [45,46]. Harp 
Lake epilimnion was sampled from March 1990 to November 
1991 for DIG concentration and ^'^G-DIG (w = 22). 

This work was initially conducted to study the origin, transport, 
cycling, and fate of organic carbon from soils to DOG and DIG in 
a soft water catchment by combining ^''G and A'''G analyses 
[28,29,42,43,45,46]. DOC turnover rates in streams, lakes, and 
wetlands are fast, <40 yr, and on the same tune scale as acid 
deposition in the area. There is extensive DOG cycling in upper 
soil layers and the source of DOG to headwater streams changes 
seasonally with water table level. Soil GO2 and ^'^G- CO2 profiles 
show characteristic patterns of root respiration, decomposition of 
SOM and diffusional gradients along with modern radiocarbon 
ages. 

The Krycklan catchment is underlain by metagreywacke 
bedrock and has forests of Norway spruce iPicea abies) and Scots 
pine (Pinus sylvestris), but deciduous trees are commonly found in 
the riparian zone of larger streams. It receives 600 mm/yr 
precipitation. Mean annual air temperature is TC. Stream 
catchment areas ranged from 0.03 km^ to 67 km^ and stream 
lengths from 0.02 km to 96.5 km. Streams were first to fourth 
order and had varying amounts of wetiand cover (0% to 40%) 
[16]. 

Stream samples were collected in June, August, and November 
2006 (« = 43). AU samples were coUected in septum-capped, screw- 
topped glass vials and analysed via headspace equilibration by gas 
chromatography or on-line GC-IRMS. Shallow groundwater was 
sampled from suction lysimeters along a transect parallel to the 
lateral flow paths towards one of the headwater streams («= 12). 
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Stream temperature was measured in the field. Typical annual pH 
range of was 3.7 to 6.3 in headwaters and 5.7 to 7.4 in 4th order 
streams. First-order-stream DOC concentrations were 5.0 mgC/L 
to 40.0 mgC/L and in 4th order streams were 5.0 mgC/L to 
15.0 mgC/L [47]. 

Research at Krycklan catchment is focused on integrating water 
quality, hydrology, and ecology in flowing waters [48]. Some 
recent work includes characterizing the loss of CO2 from stream 
surfaces by stream size and season [16,47]. Streams were always 
supersaturated in CO2, greatest in the headwaters, and negatively 
correlated with pH [16]. Owing to the importance of the gas 
exchange coefficient in controlling CO2 loss rates and its 
requirement for scaling across landscapes, it was independentiy 
measured across streams in Krycklan catchment [47]. The source 
of the excess CO2 was largely explained as respired carbon being 
exported from catchment soils [49]. As part of this work, (5"C- 
DIC values can be used to identify CO2 degassing, compared with 
the labour-intensive flux measurements, scaled flux estimates, and 
to asses the upstream CO2 loss that is demonstrated to be large 
and important [18]. 

Calculations 

We built a simple, parsimonious degassing model in Matlab 
(MathWorks, Natick MA USA) for high-DIC, low-pH, silicate 
headwater streams. Only DIG concentration, pH, ^'^C-DIC and 
temperature are required to be measured in the stream since all 
other carbonate species concentrations and ^'^C values can be 
calculated. The model partitioned total DIG into GO2, HCO3", 
and CO|^ according to pH and temperature-dependent acid 
dissociation constants and equilibrium [50,51] and kinetic [52] 
fraction factors as summarized by [34]. The calculations here are 
similar to those employed in C02sys [53], seacarb [54], and 
streamC02-DEGAS [35] but do not require as many input 
parameters, many of which are not available for our catchments. 
In small headwater streams, this simple approach may be 

Table 1. Typical model input and calculated parameters. 



appropriate to suivey a large area. Larger, more detailed models 
may be required for larger and longer streams and rivers, but it 
becomes more difficult to measure all variables required to 
ground- truth such models. 

Key assumptions to this approach are that streams are nutrient- 
poor, low pH, with low community metabolic rates relative to 
stream velocity and gas exchange, and catchments have C3 
vegetation, silicate bedrock, with high DIG shallow groundwater 
reflective of the vegetation and soils. Both Harp Lake catchment 
[21,28,29,42-46,55] and Krycklan catchment [16,24,47,56-58] 
meet these assumptions. Both have low nutrient concentrations, 
pH values, and associated metabohc rates under undisturbed forest 
canopies. Stream lengths and segments are on the order of 
hvmdreds of metres and water travel times are less than an hour. 
Forests in these catchments contain only C3 vegetation, are sub- 
boreal and boreal, and are underlain by crystalline silicate 
bedrock. Groundwaters are acidic and DIG-rich, as oudined 
below. This approach is generalizable to the plethora of small 
streams that meet these criteria across the silicate bedrock areas of 
the sub-boreal and boreal forests such as the headwaters of 
catchments on the Canadian and Fennoscandian/Baltic Shields. 
Source code is provided so that additional components can be 
added to address the complexities of carbon cycling as needed if 
the basic assumptions for field sites are not met, such as larger 
streams with longer travel times, large diel temperature variability, 
the confluence of several stream sections, measured rates of in- 
stream metabolism, known patterns of groundwater discharge, or 
the influence of carbonate-rich bedrock following portions of the 
streamG02 -DEGAS model [35], which has been used in larger 
rivers. Further, model sensitivity to variability in input parameters 
may be easily assessed with Monte Carlo approaches. 

The chemical equilibrium calculations can be easily recreated in 
any modelling language. Since ionic strength of these soft waters is 
very low (conductivity is typically 15-40 fiS), the activity 
coefficients were assumed to be approaching unity. This time- 
forward model differs from [35] in that: (a) carbonate dissolution 



Parameter 


Typical Value 


Unit 


Notes 


Citation 


Typical mode! input parameters (Fig. 2 and raw data} 


DIG 


1200 


/fmol/L 


Typical measured value 




pH 


5.5 




Typical measured value 




Temperature 


4 


°C 


Typical measured value 




^"C-DIC 


-26 




Typical measured value 




Calculated parameters based on model inputs above 


Alk 


103 


;ieq/L 


Carbonate alkalinity held constant 




CO2 


1098 


/ftmol/L 


Calculated from DIC, pH, log^„, \ogKf, 




CO2,, 


25 


^mol/L 


Calculated from logA// and pcojaim 






-10.3 


/oo 


C02*-HC03~, function of temperature 


[34, 50] 




3.3 


/oo 


HC03~~C03^~, function of temperature 


[34, 52] 


ff; 


-1.2 


%0 


CO2- CO2, function of temperature 


[34, 50] 


\ogK„ 


-6.53 




C02-HC03~, function of temperature 


[68] 


\ogKi, 


-10.6 




HC03~-C03'^~, function of temperature 


[69] 


log/t-// 


-1.18 


log(mol/l7atm) 


CO2,,,,,,- CO2, function of temperature 


[50] 




380 


ppmv 


Partial pressure of CO2 in the atmosphere 




pA'., 


14.8 




Function of temperature 


[70, 71] 
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Figure 1. DIC and CO2 degassing trajectories sKiow tliat as CO2 is lost from streams, S'^C values increase at different rates for DIC 
and CO2. 5"C-DIC is typically measured, not i5"C- CO2, but the i5''C of both DIC and CO2 during CO2 loss can be modelled using well known isotope 
fractionation factors. For this example, a typical groundwater end-member was chosen with initial DIC of 1 200 nmoUl, ij'^C-DIC of — 26%o, pH of 5.5, 
and thus carbonate alkalinity of 103 fiC- CO2, but the 5"C of both DIC and CO2 during CO2 loss can be modelled using well known isotope 
fractionation factors. For this example, a typical groundwater end-member was chosen with initial DIC of 1200 /imol/L, i5'^C-DIC of — 26%o, pH of 5.5, 
and thus carbonate alkalinity of 103 /leq/L In watersheds, variation in the groundwater end-member results in a 'family' of curves with similar 
trajectories. The grey box represents the 25th and 75th percentiles of DIC and <5'^C-DIC collected from three different first- and second-order streams 
draining into Harp Lake, Ontario, Canada. Samples were collected at the mouth of each stream from July 1989 to November 1991 (« = 111). Significant 
CO2 loss (60% to 80%) must have occurred by the sampling point in each stream in order for the box to fall where it does on the degassing curve. 
doi:1 0.1 371/journal.pone.01 01 756.g001 



was not included since our catchments are on silicate bedrock; (b) 
total alkalinity calculations differ slightly; (c) DIC concentration is 
the basis for calculations rather than total alkalinity; (d) the model 
starts with a small plug of shallow groundwater discharged to the 
surface; (e) groundwater is continuously added to the modelled 
stream if there are corresponding flow measurements; and (f) the 
iterative process to find a best fit between measured stream values 
and modelled results. 

Though there is some contribution of organic acids to total 
alkaUnity [59], loss or gain of CO2 does not change organic or 
carbonate alkalinity. Here we model alkalinity based on the 
carbonate system since in high-DIC, low-pH waters the main 
cause of pH change is CO2 loss and assume that the contribution 
of organic alkalinity to pH is relatively constant and minor. 

In practice, to calculate the fraction of CO2 lost by the sampling 
point in the stream, DIC concentration, pH, and ^''C-DIC 
measurements in the stream are required. An inverse modelling 
approach is used as the time-forward model takes initial 
constrained estimates of groundwater DIC, pH, and known 
(5'^C-DIC, and allows CO2 to degas with time. The resulting 
values of stream DIC, pH, (5''C-DIC, and ^'^C- CO2 with time 
are outputs. The model was iteratively re-run with Matlab's 
fminsearch function to reduce the sum of squared errors between 
measured stream DIC, pH, and 5''C-DIC and modelled values. 
The resulting best-fits then provide the average groundwater DIC 
concentration and the amount of CO2 that was lost upstream of 
the sampling point. In this way, an average groundwater DIC 
concentration can be determined. Additionally, the groundwater 
^'^C-DIC must be determined but this value is highly constrained 
by plant and soil ^'^C values, as above, and can be confirmed with 



a number of piezometers [28,29]. Matlab code to run the model 
forward and inversely is available from the corresponding author 
and https:/ /github.com/j|jvenky/C02-from-headwater-streams. 

At each time step, a small portion of CO2 was removed from the 
modelled aquatic system via gas exchange: 

^=fcx(C02,,-C02) (1) 

where k is the gas exchange coefficient, C02eq is the CO2 
equilibrium concentration determined with the Henry's constant 
[50], and CO2 is the solvated CO2 concentration. This necessarily 
reduced the DIC concentration since: 

dPlC ^ dC02 

DIC = CO2 -I- HCO3"- -I- CO^- (3) 

The remaining DIC was then re-apportioned to its constituent 
species by determining the pH (as H^) required to hold the 
carbonate-based alkaUnity constant [36,60]: 

H+=HC03"-l-2x CO3 " OH - - Alkalinity (4) 
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Figure 2. Groundwater and stream (a) (3'^C-DIC values and (b) DIC concentrations from Kieadwater catciiments in Sweden (left) and 
Canada (right). Groundwater and stream data are from the headwater stream network in Krycklan catchment, Sweden (K) and three first- or 
second-order streams draining into Harp Lake, Ontario, Canada (H with stream numbers; [43,44]). The boxes represent 25th and 75th percentiles, mid- 
line the median, whiskers the most extreme datum not outside 1 .5 x Inter-quartlle range (IRQ), and + the data outside 1 .5 x IQR. Number of samples Is 
Indicated In parentheses. At these two sites where we have both groundwater and stream data, there is a clear increase In ^'^C-DIC values from 
groundwater to stream sampling points along with a large decrease in DIC concentration. 
dol:1 0.1 371/journal.pone.01 01 756.g002 



The <5'^C portion of the model employed well-known temper- 
ature-dependent isotopic fractionation factors between DIC 
species, and solvated CO2 and atmospheric CO2 [34]: at 4 C, 
these are kinetic gas exchange fractionation (— 1.3%o) since the k 
for '^COj is slightly slower than for '^C02, equilibrium 
fractionation (— 1.2%o) between atmospheric COj and solvated 
CO2, CO2-HCO3" fractionation (-10.3%o), and HCO3"- 
C03^~ fractionation (-h3.3%o). At each time step, the S^^C value 
of each species was calculated. Typical values of input parameters 
and calculated parameters are summarized in Table 1. 

The effects of continuous groundwater input to the modelled 
stream that significantly changes the stream volume over the reach 
of interest can also be int lud(;d by adding an amount of water with 
known or estimated initial groundwater values (DIC, ^'^C-DIC, 
pH, alkalinity) at each time step. The amount can be determined 
by stream length and discharge measured at the catchment outlet 
or as the difference in flow between two sites. This comprised an 
additional input parameter used during best-fit modelling of the 
H4-21 data since this is a small first-order stream, the stream 
length is known, and the discharge was measured at the stream 
mouth along with DIC, ^'^C-DIC, and pH. Here, groundwater 
input was assumed to be constant down the reach. In each case, 
the CO2 lost was multiphed by the measured daily discharge and 
reported relative to the catchment area. 

Permissions 

Headwater streams in the Harp Lake catchment were studied 
under the auspices of pre-existing Ontario Ministry of the 



Environment research program. The Krycklan catchment is a 
part of the Svartberget LTER site run by the Swedish University 

of Agricultural Sciences and the area is developed for scientific 
purposes. Specific permissions for these activities were not 
required. No endangered or protected species were involved in 
either site. 

Results and Discussion 

As CO2 is initially lost from a stream, there is little associated 
change in ^"C-DIC (Fig. 1). Only after about half of the DIC has 
been lost is there an observable change of about 2%o in (5"C-DIC. 
In streams where measured (5'^C-DIC values are significandy 
different than shallow groundwater, a large amount of CO2 must 
have already been lost from the stream by the time the sample was 
collected (Fig. 1). Thus, CO2 flux measurements obtained with 
data from an individual sampling site where the 5'^C-DIC value is 
several per mUle greater than the groundwater (5'^C-DIC value 
must underestimate stream and catchment CO2 loss rates since the 
stream surface with higher CO2 concentrations and higher flux 
rates is upstream of such a sampling site. The degree of 
underestimation needs to be better quantified because these 
results suggest that net ecosystem exchange and C storage may be 
lower than assumed in boreal and sub-boreal ecozones. 

We apphed these ideas to a suite of headwater stream data from 
Krycklan catchment in Sweden and Harp Lake in Canada for 
which both groundwater and stream measurements were avail- 
able. Stream DIC concentrations were lower and ^'^C-DIC and 
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Figure 3. Typical DIC and CO2 degassing trajectories (grey band) sKiow that small initial increases in <3'^C-DIC values signify large 
losses in CO2 from groundwater to streams. The (5"C progress along this trajectory is a natural result of CO2 degassing, not other processes 
(e.g., primary production, mineral weathering, etc.). The grey band represents the combined trajectories of degassing curves with a combination of 
Initial DIC concentrations (1200 to 1800 /jmol/L), (5"C-DIC (-26%o) and pH values (4.0 to 5.5). The grey band also describes the landscape-scale 
trajectory of DIC concentrations and 5"C-DIC values from shallow groundwater discharging Into headwater streams and degassing while flowing 
downhill to small lakes. Groundwater, stream, and lake data from the two study sites are displayed as 25th and 75th percentiles (lines) with medians 
(symbols). Harp Lake data are from the epilimnion from IVIarch 1990 to November 1991 (« = 22).The line represents the degassing trajectory that was 
fit to the median of the Krycklan (open symbols) dataset. The secondary y-axis indicates the amount of DIC lost along the degassing model fit for the 
Krycklan site. 

doi:l 0.1 371/journal.pone.01 01 756.g003 



pH values greater than their respective groundwater values in both 
catchments (Fig. 2). The range of ^'^C-DIC in shallow ground- 
water, — 26%o to — 24%o is tightly constrained by the narrow 
range of ^'^C in C3 plants and soil. This is comparable to 
measured ^'^C-DOC values [28] so in-stream respiration would 
produce DIC with the same S C-DIC value as shallow 
groundwater. Surface water was allowed to lose CO2 via gas 
exchange and the chemical and isotopic equilibria were adjusted 
accordingly (see Methods). 

There are many groundwater discharge patterns possible. In 
small streams, the dilference in discharge measured at two 
locations may be the only way to estimate groundwater inputs. 
As such, it may be difficult to parameterize [43]. In a stream 
network in Alaska, high resolution CO2 and discharge data suggest 
variable and varying groundwater inputs [19] along the lengths of 
the streams. The net effect of groundwater input to small streams 
is to increase stream DIC concentration, reduce the f5'^C-DIC 
value, decrease pH, and make the stream appear closer to the 
initial shallow groundwater values (Figs 1 and 2). 

While shallow groundwaters in boreal forest catchments vary, to 
some extent in DIC concentration, (5''C-DIC, pH, and temper- 
ature, the effect of changing these variables does not alter the 
fundamental relationship between (5'^C-DIC and DIC during 
CO2 degassing. Increased initial DIC concentration, lower pH, 
and lower alkalinity cause the sigmoid-Uke relationship to 
approach its plateau (5''C-DIC value more quickly and thus 
requires a greater loss of CO2 before d C-DIC values wiU 
increase appreciably. This is demonstrated in the range of 
landscape-scale degassing trajectories (grey area in Fig. 3). This 
range is akin to a sensitivity analysis by confirming the shape of the 
DIC vs ^'^C-DIC relationship across a wide range of DIC 
concentrations and pH values. Buffered groundwaters discharging 



into headwater streams exhibit a shallower curve than acidic, 
poorly buffered waters. The degassing trajectories in DIC vs ^'^C- 
DIC space describe the CO2 loss from groundwaters to streams 
and ultimately to headwater lakes as decreasing DIC, increasing 
(5"C-DIC, and increasing pH (Fig. 3). 

Scahng-up fltjxes estimated via these methods to catchment- 
and landscape-scale requires areas or lengths of headwater 
streams, catchment areas, or groundwater discharge areas. Here, 
we avoid using stream length and surface area, parameters not 
easily obtained in headwater catchments, but see [18,20], and 
instead report stream losses relative to their catchment area. In this 
manner, the CO2 loss from aquatic surfaces is easily compared in 
the same units used to assess net ecosystem exchange and 
productivity at the catchment scale. 

One of the study catchments, a small first-order upland stream 
in a sub-boreal temperate forest (FI4-21) [61], lost a median 
15 mgC/m^/d between its source and mouth (« = 21, range 3— 
100 mgC/m^/d). This is a flux-weighted loss of 5 gC/m^/yr 
(range 1-40 gC/m^/yr) and is comparable to the annual DOC 
export (1—8 gC/m^/yr) from the catchment (Fig. 4). It is also 
around 10% of catchment net CO2 uptake [62]. Unlike DOC, this 
CO2 loss occurred from the small stream surface and was lost 
direcdy to the atmosphere before the stream outlet. Furthermore, 
this CO2 loss would not have been quantified at the mouth of the 
stream using conventional measures of CO2 and DIC fluxes [38]. 
Previous estimates of CO2 fltrx rates from the mouth of Harp Lake 
streams [2 1] are much smaller than our estimate. This highlights 
the fact that measured CO2 concentrations decline down the 
length of streams and CO2 loss is underestimated when it is 
calculated from the samphng points of lowest measured concen- 
trations. 
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Figure 4. Schematic drawing showing the size of some of the organic carbon pools, carbon fluxes, and important processes 
affecting DOC, DIC and CO2 in terrestrial catchments and aquatic surfaces. Decomposition from soil organic matter Includes microbial 
exudates. Production of DOC, decomposition, microbial uptake/sorptlon, root respiration, mineral weathering, and flushing are competing processes 
affecting the export of DOC and DIC via streams In forested catchments. Loss of CO2 directly to the atmosphere from the surfaces In the stream 
corridor can exceed the sum of DIC and DOC export to downstream lakes. Total losses of carbon lost by CO2 emissions and dissolved carbon export 
(DIC+DOC) can be Important relative to net CO2 uptake by forests from the atmosphere. All rates are per catchment area. {Figure after [42]; A. [62]; B. 
[38]; C. range for boreal and temperature forests [66]; D. [38,67]; E. this study.) 
dol:1 0.1 371/journal.pone.01 01 756.g004 



The other study catchment, a typical stream network in 
Krycklan, Sweden [18], lost upwards of 13-28 mgC/m^/d. This 
estimate is necessarily larger than the CO2 concentration based 
estimated loss of 5.0 gC/m^/yr. It is also larger than the DOC-I- 
DIC export of 4.6 gC/m^/yr and is 14-29% of estimated the net 
ecosystem exchange of 96 gC/m^/yr [18]. 

This method of using measured in-stream ^'^C-DIC values 
shows that the CO2 loss from aquatic surfaces in these boreal and 
sub-boreal catchments is large and under-estimated. The areal flux 
rates we present for Harp Lake and Krycklan catchments are 
similar to those arrived at using much more intensive sampling 
and different sets of assumptions. This supports the use of a 
parsimonious model in similar locations. The pattern of measured 
(5'^C-DIC values in small streams being much higher than 
expected shallow groundwater <5"C-DIC values, is present in 
other datasets [23,25,40] and indicates there is a large flux of CO2 
from headwater streams that has not been included in our 
continental carbon budgets. 

With a measurement at one point in a small stream, the CO2 
loss can be estimated by this ^''C approach for a certain distance 
upstream. In H4-2 1 , using a first-order-loss-rate expression QU jk, 



where U is stream velocity and k is gas transfer coefiicient, [63]), 
typical gas transfer coefficients (7-25 d^', [21]), typical measured 
mean discharge (1-20 L/s), typical cross-sectional areas (250- 
1000 cm^), the upstream distance over which the COg degassing 
occurred was 375-1375 m. 

Grab samples can be used to demonstrate the degree of 
degassing that has occurred at a given lo[:ation with a parsimo- 
nious model — the simplicity of the DIC-alkalinity-pH-(5'^C- 
DIC calculations means this idea can be easily incorporated into 
site-specific calculations. Further, stream CO2 loss, normalized to 
catchment area, is likely to be larger than CO2 loss from lake 
surfaces since lake CO2 loss is typically less than that half of the 
DOC input [5,38,64]. 

Headwater stream CO2 loss is a reduction in the net ecosystem 
productivity of boreal and sub-boreal forests. It is a globally 
important flux since boreal and sub-boreal ecozones are so large 
[1]. Here, we have shown it can be easily estimated at the 
catchment scale by combining 5"C-DIC measurements with in- 
stream and groundwater samples. Ultimately, C loss via stream 
degassing may be required to be integrated into northern 
hemispheric CO2 uptake and loss rates [2,7,18,65]. 
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